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\ Recent Hi emission observations of the Southern Galactic hemisphere have revealed 

fN| ' a mushroom- like structure extending from z = —70 to —450 pc, composed of a stem and 

a cap. Similar structures occur in three-dimensional simulations of a dynamic galactic 
J> ■ disk driven by isolated and clustered supernovae. Using these simulations, we show that 

hot gas in the Galactic disk that is not evacuated through chimneys expands into the 
■^J- , cooler gas of the thick disk, forming mushroom-shaped structures. This new class of 

objects traces buoyant flow of hot gas into the thick disk. 
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1. Introduction 

. ^ | Hi emission maps from the Canadian Galactic Plane Survey (Taylor 1999, Higgs 1999) reveal 

a structure resembling a mushroom, with a stem and a fragmenting cap (English et al. 2000). 
Assuming a distance of 3.8 kpc, the stem has a length of ~ 230 pc and the cap has a width of 
~ 175 pc. The top of the cap lies ~ 450 pc below the Galactic plane. This feature covers a total 
velocity width of ~ 24 km s _1 , with the cap redshifted from the stem by ~ 5 km s . English et 
al. (2000) estimate that the Mushroom has a kinetic energy of ~ 2 x 10 50 ergs and an Hi mass 
of at least 1.55 x 10 5 Mq. The column density contrast between the Mushroom and the ambient 
medium is no more than five, with the Mushroom having N# ~ 6 x 10 20 cm~ 2 . 

Mashchenko et al. (1999) and English et al. (2000) used thin-shell and two-dimensional hydro- 
dynamical simulations to study the origin and dynamics of the Mushroom. They modelled stellar 
winds and supernovae (SNe) in a smooth, stratified background atmosphere, and concluded that 
the observational data is consistent with a single off-plane SN having a total energy of at least 
~ 10 51 erg s -1 . In their model, the hot gas produced by the single supernova remnant (SNR) 
expanded buoyantly into a smooth gas disk. 
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In this Letter, we use the global models of Avillez (2000) to show that such mushroom-shaped 
structures can result from the buoyant rise of bubbles out of pools of hot disk gas created from 
both isolated SNRs and mergers of isolated SNRs. 

2. Supernova-Driven ISM Modelling 
2.1. Simulations 

Large-scale modelling of the interstellar gas in the disk and halo has been carried out by 
Avillez (2000) with a three-dimensional hydrodynamical code that uses adaptive mesh refinement 
(AMR). The scheme uses the piecewise parabolic method (PPM) of Collela and Woodward (1984) 
in combination with the AMR scheme of Berger & Colella (1989), combined with the subgridding 
scheme of Bell et al. (1994). 

The model includes a fixed gravitational field provided by the stars in the disk, and an ideal-gas 
equation of state. Radiative cooling is treated assuming collisional ionization equilibrium, using a 
cooling function following Kahn (1976) for 10 5 < T < 5 x 10 6 K, Dorfi (1997) for 5 x 10 6 < T < 10 7 
K, and a fit to the curves of Shapiro & Moore (1976) for the remaining temperatures. Thermal 
conduction is not included due to the numerical complexity and the lack of strong observational 
evidence for its strength in the ISM (e.g. Mac Low 1999). Its inclusion may well change our 
quantitative results, but we believe that our qualitative results on the behavior of hot gas in the 
ISM will continue to hold. 

The computational domain contains a section of the Galaxy with an area of 1 kpc 2 and vertical 
extension from -10 to 10 kpc. The innermost edge lies 8.5 kpc from the Galactic centre. The 
computational grid has a resolution of 10 pc, except in the layer between -250 and 250 pc, where 
three levels of AMR are used, yielding a finest resolution of 1.25 pc. 

Type lb, Ic, and II SNe are set up at the beginning of their Sedov phases, with radii deter- 
mined by their progenitor masses, at a rate compatible with observations (Capellaro et al. 1997), 
normalized to the volume under study. Sixty percent of the SNe are set up within associations 
(Cowie et al. 1979), and the rest are set up at random sites. 

2.2. Global Evolution of the ISM 

The disk and halo gas is initially smoothly distributed in hydrostatic equilibrium with the 
gravitational field, with scale heights given by Dickey & Lockman (1990) and Reynolds (1987). 
SNe start occurring at time zero, and the system is evolved for up to 1 Gyr. After the first 200 
Myr, the system reaches a statistically steady state that shows a dynamic equilibrium between 
upward and downward flowing gas. 
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Once disrupted by the SN explosions, the disk never returns to its initial state. Instead, 
regardless of the initial vertical distribution of the disk gas, a thin disk of cold gas forms in the 
Galactic plane, and, above and below, a thick inhomogeneous gas disk forms. The code does not 
explicitly follow ionization states, but we can trace gas with temperature T < 10 4 K and scale height 
of 180 pc, which we designate Hi, and gas with 10 4 <T< 10 5 K and scale height of 1 kpc, which 
we designate Hn. These distributions reproduce those described in Dickey & Lockman (1990) and 
Reynolds (1987), respectively. The upper parts of the thick Hn disk form the disk-halo interface, 
where a large scale fountain is set up by hot ionized gas escaping in a turbulent convective flow. 

The thick gas disk is punctured by chimneys that result from superbubbles occurring above 
the plane of the thin disk at \z\ ~ 100 pc. As they grow, they elongate along the z-direction, 
owing to the local stratification of the ISM. Chimneys in the simulation typically have widths of 
approximately 120 pc. They inject high temperature gas directly from the Galactic disk into the 
halo, breaking through the warm neutral and ionized layers that compose the thick disk. This hot 
gas then contributes to the Galactic fountain. 

3. Buoyant Outflow Prom the Galactic Plane 

About forty percent of the hot gas generated in the disk does not flow through chimneys. 
This gas still has too much energy to be held gravitationally to the stellar disk though, so it 
expands buoyantly into the thick gas disk, accelerating through the cooler gas in the thin disk, 
with characteristic velocities of some 10 km s _1 . As the hot gas rises, it entrains surrounding 
material to form a stem, often narrow, whose width does not depend on the scale height, as shown 
in Figure 1. The resulting structure resembles a mushroom composed of a cap and a stem. The 
cap drives the motion of the mushroom. During the first million years of evolution the mushroom 
consists of a rising low density cavity. The ascending gas comes from single SNRs as well as from 
pools of hot gas formed in the disk by isolated SNRs that merged together. On average, the hot 
gas produced by isolated SNe flows upward at a rate integrated over the disk of ~ 3 M & per year 
on each side of the Galactic plane, as measured at z = 140 pc (Avillez 2000). 

As the mushroom comes to equilibrium with the surrounding medium and cools down, there 
is a reversal in the column density contrast between the mushroom and the surrounding medium, 
with the mushroom cap sometimes going from a negative to a positive contrast. However, the 
positive column density enhancement is rather modest, making it very difficult to observationally 
identify mushrooms. 

Figures 2 and 3 show two mushrooms that developed at different points in the simulation. Both 
structures cooled to ambient temperatures and have column density enhancements. On average the 
ambient column density is 3 — 3.2 x 10 20 cm~ 2 , while that of the mushrooms is 3.4 — 4 x 10 20 cm~ 2 . 
The caps have the largest column density enhancement (the bluish regions inside the caps), while 
the stems of the two mushrooms have lower column densities lower than their caps (light bluish 
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regions in the stems). The mushroom in Figure 2 was created by a single SN, while the larger one 
in Figure 3 was created by three SNe. 

We can roughly measure the mass ratio rj of the cap to stem of a mushroom by estimating 
the area A and average column density N of each, and taking the mass M = AN. The mushroom 
shown in Figure 2 has a cap with A ~ 8000 pc 2 and N = 3.3 x 10 20 cm -2 ; measuring the stem 
similarly yields a mass ratio rj ~ 2.4, and a total mass of 3.7 x 10 Mq. The mushroom in Figure 3 
has a mass ratio of 4.4 and a total mass of 1.5 x 10 5 M . 

The observations of English et al. (2000) show the Mushroom to have a mass of 1.5xlO 5 M and 
a mass ratio of three (while a classical superbubble would have a ratio of 1/3). These observations 
are in rough agreement with the mushroom we show in Figure 3, suggesting that the Mushroom has 
been created by 3-4 SNRs. The greatest column density contrast between the mushrooms shown 
in Figures 2 and 3 and the ambient medium varies between 1.1 and 1.25, while English et al. (2000) 
detected a contrast of roughly two over the velocity range of the observed structure. 

4. Discussion and Comparison with Observations 

The simulations show that the disk gas is populated with strucures resembling mushrooms, 
resulting from hot gas buoyantly convecting through the thick gas disk. The dynamics of these 
structures resemble the expansion of nuclear fireballs in the Earth's atmosphere. On the Earth, as 
the fireball passes through the tropopause it halts its rise and spreads because of the turnover in 
the temperature profile. The thick gas disk has an isothermal distribution, so there is no equivalent 
to the terrestrial tropopause. However, the cap still halts when it expands to buoyant equilibrium 
with the surrounding atmosphere at a stabilization height h. 

Let us estimate this height analytically. The simulations produce a gas distribution similar to 
that of Dickey & Lockman (1990), but for simplicity in our analytical calculation we assume that the 
thick disk is isothermal with a temperature of T Q = 10 4 K and constant gravitational acceleration 
in the z— direction. Thus we approximate the thick disk with an exponential distribution given 
by p (z) = p o (0)ex.p(—z/H), where p o (0) is the density at the injection level and H = kBT /mg 
is the scale height of the medium. We choose somewhat arbitrarily a value of the gravitational 
acceleration inside the stellar disk g~5x 10~ 9 cm s~ 2 (Kuijken & Gilmore 1989), giving H ~ 75 
pc. 

The hot gas is lighter than its surroundings, so it rises, with buoyant force acting on it of 
F(z) = g[p (z) — Pi(z)], where p (z) and Pi(z) are the densities of the ambient medium and of 
the hot interior gas, respectively. As the hot gas rises, it maintains pressure equilibrium with the 
isothermal, hydrostatic surrounding gas, and continues to rise so long as a density difference is 
maintained. 

The temperature difference between the hot gas and its surroundings is expressed in terms of 
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the internal temperature gradient of the hot gas. Because the hot gas behaves adiabatically, its 
internal temperature gradient is given by 

with 

(§L-(^)t- 

where 7 = 5/3 is the adiabatic parameter and Tj(z) and Pi(z) are the temperature and pressure of 
the hot gas at height z. Noting that the hot gas remains in pressure equilibrium, 

P i (z)=P i (0) + P o (0)(e- z / H -l), (3) 

and integration of equation (1), with the help of equations (2) and (3), gives 

Tfc) =T i (0)e- z( ~ f - 1)/i ~ /H \ (4) 

Therefore, the density of the ascending hot gas is pi(z) = pi(0) exp[— z/{pjH)\. In these equations, 
z = refers to the injection level, rather than the plane of the galaxy. 

The gas stops rising when the buoyant force vanishes, and p (z) = Pi(z). This allows the 
determination of the stabilization point, whose location above the injection level is 

k = 5 - 75fflo <Si)' (5) 

Taking the scale height of the isothermal, exponential atmosphere that we assumed, H ~ 75 pc, 
hot gas with a density of pi(0) ~ 10~ 26 g cm~ 3 will stabilize at some 435 pc above the injection 
level. This is slightly above the value observed in the Mushroom and found in the simulations, 
most likely due to radiative cooling, as we now discuss. 

The adiabatic approximation for the buoyant rise only holds if the hot gas takes longer to 
radiatively cool than to rise. Let us assume that hot gas has its origin in the supernova remnant that 
generated the Mushroom. As the remnant comes into pressure equilibrium with the surrounding 
medium, it contains hot gas that eventually will mix with the surrounding medium with an increase 
in specific energy of 

Eo = * (2.5 x 10 14 erg g- 1 )P 12 p 26 1 , (6) 

where Pi 2 = Pj(0)/(10~ 12 dyn cm 2 ) and p 2 $ = Pi(0)/(10 -26 g cm -3 ) are the pressure and density 
of the hot gas. The adiabatic parameter for the mixed gas is then 

(2 \ 5 ^ 3 
-EA Pi(0)- 2 / 3 = 2.15 x 10 31 P 12P26 5/3 , (7) 
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and the time the gas takes to cool is (Kahn 1976) 

tcool = ^ = (8-3 Myr)P^ 3 p 26 W/9 , (8) 

where q = 4 x 10 32 cm 6 g -1 s~ 4 is a constant dependent on the mean mass of atoms and ions 
present in the plasma. This is a time much larger than the time needed for a supernova remnant to 
come into pressure equilibrium with the ISM, so the hot gas must expand upwards into the thick 
gas disk. 

The dynamical time for the hot gas to reach the stabilization point is td yn = {h/g) 1 / 2 ~ 
16.4 Myr for our usual parameters, roughly twice the cooling time. Thus, cooling does act during 
the buoyant rise of the hot gas, reducing the interior pressure and expansion during the rise, and 
therefore increasing the interior density and ultimately reducing the height of the stabilization point 
h to that found in the simulations. 

The description presented here applies to the mushroom-like structure described by English 
et al. (2000). We have shown that the structures found in the simulations resemble the observed 
structure in shape, size, vertical extension, and sign of column density enhancement. 

In their numerical simulations, Mashchenko et al. (1999) and English et al. (2000) attempted 
to reproduce the Mushroom and the buoyant rise of the bubble that originated it. In their toy 
models they used artificially low density, low scale-height (~ 60 pc) atmospheres and an initial SN 
remnant with a temperature of some 10 7 K. The simulations described in this Letter show that 
such a low scale height is not a requirement, and that structures similar to the Mushroom can 
occur in a fully dynamical, supernova-driven disk with a turbulent density distribution on average 
reproducing the observed distribution found by Dickey & Lockman (1990). 

The temperature of the rising gas in our model is initially only 10 6 K, and decreases further 
as a result of its adiabatic expansion. Thus, no observable X-ray emission is expected by the time 
a mushroom is fully developed with a column density higher than the surrounding medium. 

Two issues remain: what determines the column density contrast between the mushrooms and 
the surrounding medium; and how frequently does this new class of objects occur in the Galaxy? 
The column density contrast may result from enhanced cooling due to the compression of the 
mushroom cap producing higher densities in the cap. A more physical treatment of the cooling 
may lead to the higher observed column densities. 

The number of mushrooms that form is roughly correlated with the number of isolated SNe 
distributed in the stellar disk Assuming that the distribution of type Ib+c and II supernovae follow 
the distribution of the molecular gas, their rate decreases with the galactocentric radius (Ferriere 
1998), and therefore, the number of mushrooms will decrease as well at similar rate. Low density 
cavities with small sizes may be detectable as regions with column density deficits during the 
first stages of the expansion, but later these become very hard to detect because of the low column 
density contrast between them and the ISM. The simulations indicate that, even neglecting Type la 
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SNe, at least 40% of the hot gas in the thin disk, heated by isolated SNe, should eventually form 
mushrooms. 
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REER fellowship to M.-M. M. L. (AST 99-85392) and by the AMNH through a Coleman Re- 
search fellowship to M. A. Initial work was supported by ESO/FCT (Portugal) through contract 
PESO/P/PRO/1 189/97 to M. A.. 

REFERENCES 

Avillez, M. A. 2000, MNRAS, 315, 479 

Bell, J., Berger, M., Saltzman, J., & Welcome, M. 1994, SIAM J. Sci. Comp., 15, 127 
Berger, M. J., & Colella, P. 1989, J. Comp. Phys. 82, 64 

Cappellaro, E., Turatto, M., Tsvetkov, D. Yu., Bartunov, O. S., Pollas, C, Evans, R., & Hamuy, 
M. 1997, A&A, 322, 431 

Colella, P., & Woodward, P. 1984, J. Comp. Phys., 54, 174 

Cowie, L. L., Songaila, A., & York, D. 1979, ApJ, 230, 469 

Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 

Dorfi, E. A., 1997, in Proc. Saas-Fee 27, Computational Methods for Astrophysical Fluid Flow, ed. 
Steiner, O., Gautschy, A. (Heidelberg: Springer- Verlag) 263 

English, J. Taylor, A. R., Mashchenko, S. Y., Irwin, J., Basu, S., & Johnstone, D. 2000, ApJ, 533, 
L25 

Ferriere, K. M. 1998, ApJ, 503, 700 

Kahn, F. D. 1976, A&A, 50, 145 

Kuijken, K., & Gilmore, G. 1989, MNRAS, 239, 651 

Higgs, L. A. 1999, in "New Perspectives on the Interstellar Medium", ed. A. R. Taylor, T. L. 
Landecker, & G. Joncas (San Francisco: ASP), 15 

Mac Low, M.-M. 1999, Rev. Mex. Astron. Astrofis., Ser. Conf., 9, 273 

Mashchenko, S. Y., Basu, S., English, J., Taylor, A. R., k Irwin, J. A. 1999, BAAS, 194, 72.05 
Reynolds R. J. 1987, ApJ, 323, 118 



-8- 



Shapiro, P. R., & Moore, R. T. 1976, ApJ, 207, 460 

Tammann G. A., Loeffler W., k Schroeder A. 1994, ApJS, 92, 487 

Taylor, A. R. 1999, in "New Perspectives on the Interstellar Medium", ed. A. R. Taylor, T. L. 
Landecker, & G. Joncas (San Francisco: ASP), 3 



This preprint was prepared with the AAS IATjrjX macros v5.0. 



-9- 




100 200 300 400 500 I 100 200 300 400 500 

X (pc) X (pc) 




100 200 300 400 500 10 100 200 300 400 500 

X (pc) X (pc) 



Fig. 1. — Two-dimensional cuts through the three-dimensional density distribution from the simu- 
lations show buoyant expansion of hot, low-density gas into the thick gas disk. A mushroom-shaped 
structure is formed as the ascending gas expands above the 180 pc scale height of the atmosphere. 
Note the narrow stem that forms as the cap rises. 



-10- 




^400 450 500 550 600 



X (pc) 

Fig. 2. — Column density map of a mushroom located in the northern hemisphere after it reached 
the stabilization point. The color scale is in units of iV/10 20 cm -2 . This map was determined by 
integrating the column density along the y— direction of the three-dimensional data cube over a 
pathlength of 1 kpc. The maximum and minimum values in the scale are 3.4 and 2.6, respectively. 
There is a positive column density contrast between the mushroom and surrounding medium. 
However, this contrast is very low. 
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Fig. 3. — Column density map of a mushroom that developed in the northern hemisphere. The 
color scale are is in units of iV/lO 20 cm -2 . The mushroom evolved from negative to positive column 
density enhancement. The maximum column density in the mushroom is 4 x 10 20 cm -2 (in the 
cap) while the surrounding medium has a column density of some 3.2 x 10 20 cm . 



